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Abstract 

We calculate the relaxation rate of a scalar field in a plasma of other scalars and 
fermions with gauge interactions using thermal quantum field theory. It yields the 
rate of cosmic reheating and thereby determines the temperature of the "hot big 
bang" in inflationary cosmology. The total rate originates from various processes, 
including decays and inverse decays as well as Landau damping by scatterings. It in- 
volves quantum statistical effects and off-shell transport. Its temperature dependence 
can be highly nontrivial, making it impossible to express the reheating temperature 
in terms of the model parameters in a simple way. We pay special attention to the 
temperature dependence of the phase space due to the modified dispersion relations 
in the plasma. We find that it can have a drastic effect on the efficiency of perturba- 
tive reheating, which depends on the way particles in the primordial plasma interact. 
For some interactions thermal masses can effectively close the phase space for the 
dominant dissipative processes and thereby impose an upper bound on the reheat- 
ing temperature. In other cases they open up new channels of dissipation, hence 
increase the reheating temperature. At high temperatures we find that the universe 
can even be heated through couplings to fermions, which are often assumed to be 
negligible due to Pauli-blocking. These effects may also be relevant for baryogenesis, 
dark matter production, the fate of moduli and in scenarios of warm inflation. 
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1 The big bang temperature 



Many properties of the cosmos we observe today can be understood as the result of pro- 
cesses that occurred during the early, high temperature phase of its history [1]. This makes 
the maximal temperature in the early universe a crucial parameter in physical cosmology 
[2] . An important probe of the thermal history of the universe is given by the abundance 
of thermal relics. These include the light elements that we observe in the intergalactic 
medium. They were created in thermonuclear reactions in the primordial plasma, known 
as big bang nucleosynthesis (BBN). The good agreement between theoretical BBN calcu- 
lations and observation allows to conclude that the universe in the past had a temperature 
Tbbn of at least a few MeV. 

Observationally, only little is known about the thermal history of the universe beyond 
that. However, there are good reasons to believe that it has been exposed to much higher 
temperatures. If the observed dark matter (DM) is a thermally produced weakly inter- 
acting massive particle, then temperatures not too much below that particle's mass are 
required for efficient thermal production Also other DM scenarios often require tempera- 
tures T ^> Tbbn, for instance sterile neutrinos [3, 4] (see [5] for a recent review) or axions 
(see [6] for reviews). The observed baryon asymmetry of the universe (see [7] for a recent 
review) provides another hint for large temperatures. If the baryon number violation in the 
Standard Model of particle physics (SM) is the relevant source of the baryon asymmetry 
[8] ( "electroweak baryogenesis" ) , the universe must have been exposed to temperatures 
T > T EW ~ 140 GeV. Many models of baryogenesis require a temperature that is even 
much higher than that, including GUT-baryogenesis, thermal leptogenesis [9], baryogen- 
esis from neutrino oscillations [3, 10] or many variants of Affleck-Dine baryogenesis [11]. 
Hence, any upper bound on the temperature in the universe allows to rule out or constrain 
models of particle physics beyond the standard model (SM). A well-known example is the 
conflict between supersymmetry and thermal leptogenesis due to the gravitino problem 
[12]. A large temperature can also be related to the fate of moduli [13] and the be of great 
importance for the possible decompaticfication of extra dimensions [14]. 

1.1 The origin of the hot big bang 

It is widely believed that the hot big bang initial conditions of the radiation dominated 
universe were created during reheating after a period of accelerated cosmic expansion (cos- 
mic inflation) [15, 16]. Inflation made the universe homogeneous and flat, while creating 
the seeds for later structure formation from quantum fluctuations [17]. The accelerated 
expansion can easily be explained if one assumes that the energy density of the universe 
was dominated by the potential energy of a scalar field 0, the inflaton 1 , which leads to a 
negative equation of state while is slowly rolling towards the potential minimum. At the 
end of the inflationary epoch, the universe was cold and empty, with all energy stored in 
the zero mode of <fi. All other particles were then created and the universe was (re)heated 

■'■For simplicity, we consider a single field model of inflation. However, the considerations in section 2 
apply to any model of inflation that includes perturbative reheating. 
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[18, 19, 21, 22] by the dissipation during coherent oscillations of <fi around its potential 
minimum. The early stage of reheating, sometimes referred to as preheating, most likely 
included complicated nonequilibrium processes such as particle production [19-22], tur- 
bulence [23] and non-thermal fixed points [24] or instabilities [25, 26], including explosive 
particle production by parametric resonance [19, 22, 27, 28]. Some of these effects have 
recently been discussed in [29]. However, the late phase of reheating, which determined 
the temperature at the onset of the radiation dominated epoch, in many scenarios can be 
treated perturbatively. Thus, the hot big bang initial conditions were set by perturbative 
processes involving individual inflaton quanta. In addition, the efficiency of perturbative 
processes is also of great interest in the context of preheating as they may delay the para- 
metric resonance [30]. 

Sometimes the reheating temperature defined by assuming that a) the inflaton converts 
its entire potential energy into radiation instantaneously at the end of inflation and b) this 
radiation comes into thermal equilibrium instantaneously. In reality the inflaton dissipates 
its energy at a finite rate T, and the reheating process takes place over a time interval ~ 1 /Y. 
As the universe expands during this period and the radiation cools down, the assumption 
of instantaneous conversion leads to an overestimate of the maximal temperature in the 
universe. The temperature dependence of Y is crucial to understand the time evolution of 
the temperature. Indeed, the temperature decreases during much of the reheating process, 
and the temperature at the onset of the radiation dominated era may be considerably 
lower than the highest temperature the universe has been exposed to. The reheating 
temperature should be understood as the highest temperature in the radiation dominated 
universe [1, 31, 32]. 

The presence of the primordial plasma affects the rate Y in various ways: Quantum 
statistical factors can enhance or suppress individual channels due to Bose enhancement or 
Pauli blocking. The modified dispersion relations of particles in the hot plasma make the 
phase space temperature dependent and new channels of dissipation can open up in addition 
to the decay of 0-quanta, such as Landau damping. As we will discuss in the following, 
this may greatly enhance or suppress Y and have a strong effect on both, the maximal 
temperature the universe has been exposed to and the temperature at the beginning of the 
radiation dominated epoch. 

1.2 The role of "thermal masses" 

The inflaton mass m is thought to be much larger than the masses of known particles, 
and the SM is in the symmetric phase for T > Tew- in vacuum the only perturbative 
mechanism by which dissipates its energy into other degrees of freedom is the decay of 
heavy <p quanta into lighter states. The particles found in the universe today are the end 
products of the subsequent decay chain. However, once <j) has lost a significant fraction 
of its energy, the universe is filled with a medium formed by the hot plasma of its decay 
products. For the moment we collectively refer to these as Afj, without specification of the 
spin, charges or other quantum numbers of these fields. The Xi represent all degrees of 
freedom other than (j), including the SM fields. The inflaton relaxation time, characteristic 
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for the duration of reheating, is given by r = 1/T, where T is the <fi relaxation rate in the 
medium. The coupling of <p to other fields should be very small in order not to spoil the 
flatness of its potential required for slow roll inflation [33]. This implies that r is much 
longer than the typical time scale associated with the relaxation to thermal equilibrium 
for any other field in the plasma. 2 This suggests that it is reasonable to approximate 
perturbative reheating as an adiabatic process, during which the Xi are close to kinetic 
equilibrium at any time 3 . 

It is well-known that in quantum field theory the spectrum of propagating states 
( "quasiparticles" ) and their dispersion relations are modified in a medium due to screening 
and collective phenomena. In the following we assume for simplicity that the Xi are in ki- 
netic equilibrium and can be characterised by an effective temperature T 4 . Then thermal 
corrections to the dispersion relations typically scale as ~ a{T . Here «j are dimensionless 
numbers that characterise the strength of a typical interaction amongst the Xi. They can 
e.g. represent gauge or Yukawa coupling constants. In the simplest case, the system can be 
pictured as a gas of quasiparticles with momentum independent thermal masses Mi ~ a{I '. 

While 4> should have very weak interactions, the at T ^> Tew can be of order one. 
Thus, the thermal masses of the decay products Xi grow much faster with temperature 
than the effective <p mass. In [34] it has been concluded that due to this effect, the inflaton 
decay becomes kinematically forbidden at some critical temperature T c , when the sum of 
the effective masses of the decay products exceeds the inflaton mass, 



As a result, the universe could never be heated up to temperatures larger than T c . Due to 
the very small couplings of 0, thermal corrections to the <fi dispersion relations are tiny and 
the LHS of (1) can essentially be identified with the vacuum mass m. 5 Thus, T c , defined 
by (1), in good approximation does not depend on the inflaton coupling. If T c imposes an 
upper bound on the reheating temperature, this would have far reaching consequences as 
it would allow to constrain the temperature in the early universe without knowledge of the 
inflaton couplings. 

In [37] it was shown for a simple toy model that V is non-vanishing in the limit T T c . 
However, the analysis does not take into account the fact that V may be vanishing for 
intermediate temperatures T ~ T c ; in this case T c would pose an upper bound on T, as 

2 We do not make any other assumptions about the strength of the inflaton coupling. It may indeed be 
a non-trivial requirement that it is large enough to allow for efficient reheating and small enough not to 
spoil the flatness of the inflaton potential. We do not address this issue here. Wc express all contributions 
to r in units of their zero temperature value, in which case its dependence on the inflaton coupling cancels 
at leading order. 

3 It is well-known that the thermalisation process in bosonic toy models can be prolonged because most 
energy is stored in the infrared modes during preheating. However, in more realistic models involving 
fermions equilibration can occur considerably faster [35] . 

4 We make this assumption for computational simplicity; the arguments brought forward in section 2 
remain valid even if the background is far from equilibrium [36] . 

5 The generation of an effective mass in the radiation dominated era has e.g. been studied in [39]. 
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the universe could never pass this "forbidden regime" of temperatures. We indeed find 
that in the toy model used in [37] this is the case because T is strongly suppressed at 
T ~ T c . The problem was studied in more detail in [38], where a physical interpretation 
in terms of scatterings was given. In the present work we extend that study and perform 
a detailed analysis of the kinematics in a hot plasma and its effect on T. We show that 
the suppression of T at T c due to thermal masses only occurs in special cases; in general 
they do not impose an upper limit on the reheating temperature. We also find that at high 
temperatures the dissipation into fermions can give a significant contribution to L. 

In section 2 we argue on general grounds that the kinematics of reheating can often not 
be described correctly by simply replacing vacuum masses of particles by thermal masses. 
In sections 4-7 we illustrate this by calculating the inflaton relaxation rate T in a model of 
chaotic inflation for different interactions. We first study the dissipation into bosons. We 
chose three different examples to represent cases where thermal masses block the inflaton 
decay (section 4), do not block the decay (section 5) or Landau damping drives reheating 
(section 6). In section 7 we study the inflaton coupling to fermions with gauge interactions. 
At low temperatures the decay into fermions is suppressed due to Pauli blocking before 
thermal masses become relevant, but at large temperatures (f> can efficiently dump its energy 
into fermions via Landau damping. In section 8 we discuss our results and conclude. The 
reader who is not interested in technical details may skip sections 4-7 and read sections 1, 
2 and 8 short letter. 

2 Summary: The kinematics of cosmic reheating 

Due to the medium corrections to dispersion relations the phase space in a hot plasma is 
temperature dependent. This can have a strong effect on the efficiency of reheating. When 
the primordial plasma is modelled as a gas of quasiparticles, it is often assumed that their 
kinematics can be described by simply replacing vacuum masses of particles by momentum 
independent thermal masses ~ a{I '. In the context of reheating, such approximation may 
fail even at the qualitative level, as the following effects need to be taken into account. 

• Particle production - The inflaton field ((f)) oscillates with a frequency ~ m, which 
is higher than the mass of most other particles. In such a rapidly time-varying 
background it is not possible to define an appropriate vacuum state and particles 
[20]. In such a background dissipation is not due to perturbative processes, but to 
particle production from the background. In this regime the coupling to the (<f))- 
background can give rise to the explosive particle production when it dominates 
during preheating, and arguments based on single particle kinematics certainly fail. 
A quantitative treatment of this preheating process goes beyond the scope of this 
work. 

• Quasiparticle description - When the produced particles start making up some 
fraction of the universe's energy content, they form a thermal background, to which 
we assign an effective temperature T. The existence of well-defined quasiparticles in 
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this plasma is related to the requirement that all dressed spectral densities (to be de- 
fined later) as a function of energy u feature narrow peaks at energies u = Qj, which 
define the dispersion relations of quasiparticles. These receive three contributions: 
intrinsic (vacuum) masses rrii, a time-dependent mass due to the coupling to the os- 
cillating mean field ((f)) and thermal masses due to forward scatterings in the plasma. 
A (quasi)particle description is appropriate if the typical energy a; of a A'j-mode does 
not change considerably during a scattering process [36]. This is the case if either 
the inflaton's oscillation frequency is small compared to u or the contribution to the 
particles' effective masses from the slowly changing thermal background dominates 
over that from the rapidly oscillating (0) -condensate. This is possible even if the en- 
ergy stored in (0) exceeds the energy of the primordial plasma because the inflaton 
couples weakly: If g is the inflaton coupling and «j ^> g a typical coupling within the 
plasma, then thermal contributions ~ a{T can dominate over g(<p) even if (</>} 3> T. 
In this case reheating can be described in terms of perturbative processes involving 
individual quasiparticles and phase space arguments apply. This is the regime we 
focus on in this work. Since the (^-amplitude decreases while T increases in the 
course of reheating, it becomes an increasingly good approximation at later times. 

• Landau damping - In the regime where a quasiparticle description is appropriate, 
large thermal masses may kinematically block the inflaton decay. However, decay is 
not the only channel by which (p can dissipate energy into the plasma. 0-quanta can 
engage in scatterings with particles from the bath, leading to a dissipation mechanism 
similar to Landau damping [53]. At low temperatures, this mechanism is inefficient 
because the concentration of scattering partners is too low. At high temperatures 
it becomes the dominant channel of dissipation because processes involving several 
interactions in a given time unit become equally likely as those involving one or few 
when the density is high. Landau damping is the main reheating mechanism in the 
model discussed in section 6, but also appears in sections 4, 5 and 7. 

• Quasiparticle dispersion relations - The dispersion relations for quasiparticles 
in most cases cannot be obtained by simply replacing bare masses with momentum 
independent thermal masses. In general, the quasiparticle energy f2j is a compli- 
cated function of momentum p. In spite of that, many transport phenomena can 
effectively be described by thermal masses Mj ~ a{T . The reason is that for hard 
momenta p ~ T the dispersion relations usually asymptotically approach a form 
Vti ~ (p 2 + M 2 ) 1 / 2 . Since most particles in a plasma have momenta ~ T, a de- 
scription in terms of momentum independent thermal masses is appropriate for most 
microscopic processes. However, when the inflaton decays, the decay products have 
typical momenta < m/2 (and not ~ T). For m < T these are "soft" compared 
to typical energies ~ T in the plasma. The dispersion f2j relations for soft modes 
strongly depend on the details of the interactions; they can be complicated functions 
of momentum (see sections 7 and 5). It is not necessarily true that these modes 
become "thermally heavy" and kinematically block the inflaton decay, see section 5. 
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• Collective excitations - It is well-known [52] that for soft momenta there can be 
collective excitations in a plasma (holes/plasminos, plasmons...) that behave like 
quasiparticles in addition to the screened versions of particles known from vacuum 
[68, 69], see [70, 72] for a recent discussion. Decays into these and scatterings with 
them may give a contribution to T [54, 55], as we see in section 7. 

• Off-shell transport - In quantum mechanics processes involving virtual (off-shell) 
particles contribute to transition amplitudes. These contribute to the energy transfer 
between the inflaton and other fields even if some of the intermediate states are too 
heavy to be produced as real (quasi)particles. In our setup the relevance of off-shell 
processes can be parametrized in terms of the thermal X r widths IV Though in 
principle suppressed by I^/f^, the off-shell transport can give a sizeable contribution 
to r due to collinear amplification or large occupation numbers [37, 38, 40, 48, 50]. 

These considerations are very general. In the following sections we illustrate them in a 
simple model of chaotic inflation. The reader who is not interested in technical details may 
skip them and read sections 1, 2 and 8 as a short letter. 

3 Elements of thermal quantum field theory 

The crucial quantity that determines how efficiently <fi can reheat the universe is the relax- 
ation rate V. It determines how much energy dissipates into the primordial plasma per 
unit time. The "in-out formalism" and S-matrix, commonly used in quantum field theory 
in vacuum, do not provide an appropriate tool to describe nonequilibrium phenomena in a 
dense plasma. The reason is that a nonequilibrium process is an initial value problem, in 
which the final state is not known a priori and memory effects can be important. Further- 
more, the definition of asymptotic states is ambiguous in the omnipresent plasma when 
the density is large enough that particles always feel the presence of their neighbours. A 
consistent treatment is possible in the framework of nonequilibrium quantum field theory, 
where all properties of the system can be expressed in terms of correlation functions of 
quantum fields. 6 The approach we use in the following is known as Schwinger-Keldysh 
formalism [58], but sometimes also referred to as "closed time path" or "in- in formalism". 
We use the notation of [38, 51, 57]. 

In the Schwinger-Keldysh formalism, the gain and loss rates r< and r> for the 0-mode 
q are related to self energies n < (xi,x 2 ) and n > (xi,x 2 ). Inflaton couplings that are linear 
in <p the interaction terms can be expressed as <f)0[Xi], where 0[Xi] is represents operators 
that are composed of fields X^ other than <fi. In the following sections we will use the model 

6 For a review see [56]. 
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Lagrangian 



i=i ^ 

2 , 

-^70X1X2 - Yc/>m - + ( 2 ) 

The Xi are scalar and \l/ is a fermionic field to which <fi couples. Lx contains all other degrees 
of freedom (including the SM fields) and their couplings to the Xi an d In this section we 
do not specify the interactions in Lx] as before, we symbolically characterise their strength 
by dimensionless numbers a,. It is assumed that m > m^m and g/m,Y,hi <C a; t . For (2) 
the operator 0\X^\ is identified with 

2 , 

G[x] = g X i(x)x2(x) + YV(x)*(x) + ^jxKx)- (3) 

i=l 

To leading order in the tiny inflaton couplings the self energies can be calculated as 

n > (x 1 ,x 2 ) = (0(x 1 )0{x 2 )), U<( Xl ,x 2 ) = (0(x 2 )0{x 1 )), (4) 

The average (...) is defined in the usual way as (^4) = Tr(gA), where g is the density 
matrix of the thermodynamic ensemble under consideration. It includes the usual quantum 
average as well as a statistical average over initial conditions. It is convenient to define the 
spectral self energy, 

Yl~(xi,x 2 ) = n > (xi,x 2 ) - n < (xi,x 2 ), (5) 

which can be related to the usual retarded self energy by H R (xx, x 2 ) = 9(ti — t 2 )n~(xi, x 2 ). 
It is sometimes also referred to as dissipative part of the self energy, as it determines the 
total relaxation rate for 0, including gain- and loss terms [53], cf. (15). 

In the following we will for simplicity always assume that all fields other than are in 
thermal equilibrium with temperature T (they relax on time scales l/T^ 1/r, while we 
are interested in the time scale 1/r, which is characteristic for the duration of reheating 7 . 
Due to the homogeneity and isotropy of the universe, IT - depends only on the relative 
spacial coordinate |x x — x 2 |. To leading order in the tiny couplings g, hi and Y, the self 
energy Il~ is computed from diagrams that have no internal 0-lines. 8 Since all other fields 
are in equilibrium, Il~ to this order is invariant under translations in time and depends 
only on the relative coordinate X\ — x 2 . In principle it is straightforward to generalise the 

7 In reality thermalisation is a highly non-trivial process. However, most of the following can be gener- 
alised in a straightforward way as long as T <g; m [36] . 

8 Diagrams involving the mean field (or one point function) (0), as shown in figure 1, are of higher order 
in the tiny couplings g, hi and Y . This suppression is, however, compensated when the amplitude of the 
(</>)-oscillations is bigger than ~ onT 2 jg. We exclude this case, which may be realised during the early 
stage of reheating, as we are interested in the perturbativc regime. 
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Figure 1: An example for a contribution to II that involves the mean Geld; crosses represent 
the mean held (<p), lines represent propagators for <fi (solid), xi (dashed) and \i (dotted). 



following considerations to a time dependent background [36], but the computation of T 
in practice becomes much more difficult. Thus, we can define the Fourier transforms 

nSC) - /✓<*, - *) ^e-^-xnn*, - „). (6) 

Another quantity that we will need in the following is the spectral function A~( X i,x 2 ), 

A~(x 1 ,x 2 ) = i ((0(xi)0(x 2 )) - (0(x 2 )0(xi))) . (7) 

When n~ is translation invariant, the same applies to A - [57], and we can define the 
spectral density p q (w) as 

p q (w) = -i J d\x x - x 2 ) e laj(tl -* 2) e- iq{xi ~ X2) A-(xi - x 2 ). (8) 

It can be expressed as [57] 

( ] = -jlmJlgH + 2ue 

P ^ U) (cu 2 Ren«(cu)) 2 + (Imn*(w) + cue) 2 ' 1 ' 

The pole structure of p q (cj) in the complex u plane determines the spectrum of propagating 
states in the plasma. Let fi q be a pole of Pq(u) with f2 q = Ref2 q and T q = 2Imr2 q . Both 
are temperature dependent because Ilq (u) depends on T. In weakly coupled theories one 
usually observes the hierarchy 

r q < n q . (io) 

It that case p q (w) features peaks of width ~ T q at energies u ~ ±^ q which can be 
interpreted as quasiparticle-reson&nces. They may correspond to screened particles or 
collective excitations of the plasma 9 . It is clear from (9) that, if (10) is fulfilled, the 
quasiparticle dispersion relation (or "mass shell") fi q is essentially fixed by Relived) via 
the condition 

u 2 - q 2 - m 2 - RenJ(cu) = 0, (11) 



9 Some authors use the expressions "quasiparticle approximation" and "narrow width approximation" 
to refer to the limit r q /f2 q — > 0, which corresponds to a spectral density (20). We use these terms in 
a broad sense to refer to the situation of a small, but finite width described by (12), corresponding to 
unstable, but long-lived quasiparticles. 
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while Imllq(u;) gives the thermal width. The dispersion relation f2 q in general has a 
complicated q-dependence. Only in special cases it can be parametrized by momentum 
independent "thermal masses" 10 . RelT^u;) contains a zero temperature divergence that 
can be absorbed into the physical mass in the same way as in vacuum [57, 59]. In the 
following we interpret m as physical mass, with Rell^ (u) finite. The physical properties of 
quasiparticles can be read off a Breit-Wigner fit to (9) near the pole, which is parametrized 



by T q and f2 q . 



BW/ \ 
Pq M 



2Z, 



%) 2 



(12) 



cont , 



is the continuous part of p q (w). To obtain the correct residue, we introduced the 



parameter 



Using the relation 



1 - 



2a 



n; 



<9ReITj(w) 
du 



2iImnJ(o;) 



(13) 



(14) 



we can express the thermal width as 



Zr 



2fi q " q 



(15) 



For the (ft width we can neglect the factor Z^ because of the small inflaton coupling. For 
the other fields Xi it should be kept, cf. (30). T q yields the total relaxation rate (or 
damping rate) for [57], which plays the role of a production rate when the occupation of 
the mode q is below its equilibrium value and that of a dissipation rate in case it is above 
its equilibrium value, as in the case of reheating. T q can be decomposed into gain- and 
loss rates r q = r> — T^, with T| = ill|(f2 q )/(2f2 q ). The relaxation for the coherent zero 
mode (eft) is given by T = r , see e.g. [57, 59, 60, 62, 63]. We assume that due to the very 
weak inflaton coupling, thermal modifications of the <j) dispersion relation are negligible. 
Then there exists only one type of eft quasiparticle with f2 q ~ w q = (q 2 + m 2 ) 1 / 2 . On the 
other hand, corrections to the thermal 0-width T q are important. 

In the models discussed in the following sections <fi couples to other scalars Xi an d 
fermions \1/ described by (2). The spectral densities for Xi an d are defined in analogy 
to (7) and (8), 



PipiPo 

>p(Po))q 



d\ Xl - X2 ) e *o(ti-*0 e -ip(*i-«) ((xiixjxiixz)) - {Xi{x 2 )xi{xx))) (16) 
d\ Xl _a; 2 ) e w(ti-t 2 ) e -ip(xi-x 2 ) (^ Q ( Xl )^( X2 )) + (^( X2 )M/ Q ( Xl ))) , 

(17) 



10 Such a parametrization can be a good approximation if only a certain range of momentum modes is 
relevant. For very soft modes one may use the plasma frequency, i.e. the solution to (11) for q = 0, for 
hard modes q ~ T one often uses the asymptotic mass. 
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Figure 2: Relevant Fcynman diagrams, solid lines represent 4>, dashed lines xi, dotted lines 
Xi, solid lines with arrow ty. 



where a and p are spinor indices, which we suppress in the following. The different sign 
in the fermionic spectral density originates from the Grassmann nature of the fields. Self 
energies II~ and S~ for the x% an d ^ can be defined in analogy to U^, but of course involve 
all interactions of the field under consideration. 

In scalar theories at high temperature, there is usually only one type of quasiparticle 
that is relevant for transport [52]. This justifies the (one) pole approximation 

pole/ \ _ ~ ^Ap(Po) , , 

In the simplest case the x% dispersion relations can be parametrised as fl ip = (p 2 + 
M^T)) 1 / 2 , with momentum independent thermal masses Mj(T). We will also use the 
Breit-Wigner approximation 

and its "zero width limit" , 

p°,(Po) = 27rZ ip sign(p )5(p 2 - fi 2 p ) + p™\p ). (20) 

Here f2j P = Ref2jp and T ip = 2Imf2j P are the quasiparticle energy and thermal width for Xi 
excitations with spacial momentum p. The continuous contribution due to multiparticle 
states is often suppressed in weakly coupled theories, see [59, 63] for some discussion. 
Therefore it is often neglected. The fermionic spectral density p p {po) is more complicated, 
we discuss it in section 7. 



4 Example I: blocking by thermal masses 

We first discuss the dissipation into bosons. We consider a situation where the interactions 
of can effectively be described by the following model 

i=l \ ' 

(21) 
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Here (j) is the inflaton, Xi are t wo other scalars with rrii <C m and £bath represents all other 
fields in the primordial plasma. Using the real time formalism of thermal field theory (see 
e.g. [52]), the contribution to T from the diagram shown in figure 2a) is obtained from 11 




1 (1 + /b(Po) + /b(w - Po)) Pi P (Po)p2q-p(w - p ). (22) 



Here /b(o;) = (e^/ T — 1) 1 is the Bose-Einstein distribution. We will use different approx- 
imations to the spectral densities p». 

4.1 Zero width approximation for quasiparticles 

The structure of (22) is general for leading order diagrams. The product of spectral den- 
sities for the internal lines sorts out the kinematics of contributing processes, while the 
distribution functions account for the stochastic weight in the ensemble. When plugging 
in (20) for the spectral density, (22) is only non-zero when both of the particles in the 
loop can be on-shell simultaneously. This is the case when the submanifolds defined by 
the 5-functions in (20) intersect somewhere in the integration volume. It is in full analogy 
to the optical theorem and cutting rules at zero temperature [61]. In vacuum, the total 
decay rate of a particle can be computed by cutting the self energy in all possible ways 
and requiring that all cut lines must be on-shell. The different cuts are then interpreted 
as physical processes contributing to the decay rate, with the cut propagators as outgoing 
particles (decay products). The difference at finite temperature is that the cut propaga- 
tors can also act as incoming particles, cf. figures 3 and 4. The physical reason is that 
the medium provides scattering partners for <fi. In the Schwinger-Keldysh formalism, these 
processes are automatically taken into account in the calculation of n~. Integration of (22) 
with (20) for q = yields the known expression 



with Qi = Q p i and fi = /e(fij). We have neglected the continuum part in (20). The 
different 5-functions in (23) arise from different cuts through the self energy shown in 
figure 2a) [53]. They are displayed in figure 4. The term proportional to 5(u — fii — f2 2 ) in 
the first line describes decays and inverse decays 4> -H- X1X2, see figure 4a). The two terms 
in the prefactor can be identified with the thermodynamic probabilities (j\ + l)(/2 + 1) for 
a decay (the /j are due to Bose enhancement as the Xi are bosons) and /1/2 for an inverse 
decay; they obey the detailed balance relation. The 5(u + Q± + f^) term comes with the 




x 



((/1 + 1) (/a + 1) - /1/2) (S(u - fii - n 2 ) - 6(cj + Q x + ft 3 )) 

+ + 1)/ 2 - (/ 2 + (6(u - fii + n 2 ) - 5(uj + fix - n : 



(23) 
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Sec e.g. [38; 
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Figure 3: The self-energy II = — 2«ImII can be calculated from the sum of all cuts at a 
given order. At finite temperature, cut propagators can correspond to incoming or outgoing 
lines [61]. For real scalar fields, a cut propagator that acts as incoming particle carries a 
factor Jb (representing the abundance of scattering partners in the plasma), one that acts 
as outgoing particle has a factor 1 + /b (leading to Bosc enhancement) . The cut shown in 
this figure, for instance, includes 2 O 3 scatterings. Solid lines represent <p, dashed lines xi, 
dotted lines xi ■ 




Figure 4: Different cuts through the leading order contribution to U from the trilincar 
interaction; solid lines represent <f>, dashed lines xi, dotted lines X2- 

same prefactor, but only contributes for uj < 0; it represents the generation of 4>XiX2 from 
the vacuum and the inverse process, which give no on-shell contribution. The terms in the 
second line can be interpreted correspondingly: they represent processes X2 <PXi an d 
Xi -H- (pX2, see figure 4c) and d). 

For the quartic Xi~ se ^ interaction in (21) the Xi-dispersion relations in the plasma can 
at leading order in Aj be parametrized as fli = (p 2 + Ml)^ l ^ 2 \ with momentum independent 
effective thermal masses that are generated by diagram shown in figure 5a) 12 

Here we have suppressed the explicit temperature dependence of Mj(T) 13 . We assumed 

12 The thermal mass correction here is ~ A 1 / 2 T rather than ~ AT because it arises from a local (tadpole) 
diagram. 

13 In the following we always use capital letters for effective masses and small letters for vacuum masses, 
often without making the dependence on T explicit. 
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3) 



b) 



Figure 5: Diagrams contributing to the Xi self energy to order Xf. 



that the zero temperature part of the tadpole diagram has already been absorbed into 
rrii by an appropriate renormalization procedure 14 . The momentum independence of the 
diagram in in figure 5a) also implies Z$ = 1 at leading order in Aj. 
Then the integral (23) can be evaluated analytically [57, 59] as 



with 



8lY W 2 + Wi 



(25) 



(1 + f B ( Wl ) + f B (w 2 ))6(u - Mi - M 2 )6(u) (26) 



[9], 



Mi 



57T W 2 — W\ 



(/bW - /bM)0(-w - Mi + M 2 )9{u) 



where we have used the abbreviations 



Wi 



J 1 + M\ - Ml 
2cu 



w 2 



J 1 - Ml + Ml 
2qj 



(27) 



(28) 



The terms have clear physical interpretations. T)q\uj) is the total contribution from the 
decays and inverse decays (p -H- XiXz- The ^-function shows that such processes are only 
allowed if M ^ Mi + M 2 , as suggested by on-shell quasiparticle kinematics. Sq(ui) is the 
contribution from processes -H- Xji which give rise to Landau damping. It comes with 
a ^-function that requires |Mi — M 2 \ ^ M. In between there is a kinematic region in which 
neither of these processes is kinematically allowed. 

Let us summarize this behaviour on more general grounds. Let us assume that there 
is an interaction that gives thermal masses Mi ~ a{T to the decay products, where ai is 
a dimensionless small parameter that we use symbolically to indicate scales in a weakly 
coupled system. For a particular interaction it has to be related to the coupling constant. 
For instance, to parametrize the thermal mass (24) as c^T, we have to identify = 
(Aj/24) 1 / 2 , see (24). 15 . Since thermal corrections to the inflaton dispersion relation are 



14 The dispersive part of the self energy can only be absorbed into m, at a particular T and appears 
explicitly for all other temperatures. It is common to impose renormalization conditions at T = 0. Then 
rrii has the meaning of a physical mass in vacuum. The T-dependent part of the dispersive self energy is 
finite and no counter terms in addition to the T = case are needed [52] . 

15 Thc unusual square-root scaling is due to the tadpole diagram in figure 5. For fermions with hard 
momenta one has to identify on = a/2, see (62). 
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negligible and its width is small, we can set M ~ m. The Xi vacuum masses are much 
smaller than M ~ m. However, due to CKj ^> g/m, they grow much faster with T and we 
can define critical temperatures 

m 2 = (M 1 (T C ) + M 2 (T C )) 2 , m 2 = {M 1 {f c ) - M 2 (f c )) 2 . (29) 

We can estimate the critical temperatures as T c ~ — ^ — and T c ~ , — — — r . For T < T c , 

(p can reheat the plasma by decay, for T > T c by Landau damping. If it can be justified 
to use (25) to compute the total relaxation rate T, T c indeed poses an upper bound on 
the temperature because the primordial plasma can never reach temperatures T > T c , 
where Landau damping would be efficient. However, with the Xi self-interactions (21), the 
trilinear coupling gcpXiX2 allows for scatterings (f>Xi ^ XjXjXj? see figure 3. These require 
the exchange of an intermediate Xj quantum 16 . In the region T c < T < T c , this intermediate 
state must be off-shell. In addition, there is a suppression by the coupling constant. In 
spite of that they can be relevant at high temperature because both of these suppressions 
may be compensated by the large occupation numbers in the plasma. Whether or not the 
approximation (25) is justified has to be decided on a case by case basis. 

4.2 Inclusion of thermal widths and off-shell transport 

In the Schwinger-Keldysh formalism, all contributions are taken into account consistently 
when the 0-self energy U~(u) is calculated at higher order in the couplings aj. A consistent 
computation of n~ (cj) to a given order in all couplings automatically includes all processes 
that contribute to T at that order. This may require resummation [64] to fix the so-called 
breakdown of perturbation theory 17 at high temperature [52]. The need for resummation is 
usually related to infrared or collinear divergences. Infrared problems at finite temperature 
mathematically are related to the divergent behaviour of the Bose-Einstein distribution 
function fg for small arguments. Physically, they can be understood easily: At high 
temperatures, the density of potential scattering partners in the plasma becomes so large 
that the probability for multiple scatterings within some time unit is of the same order of 
magnitude as that for a single interaction. Collinear enhancement in thermal field theory 
[50, 75] can occur when the decaying particle is relativistic and the decay products move 
parallel for a long time in the rest frame of the bath, giving them enough time for multiple 
exchange of force carriers. 

For the Lagrangian (21) the full relaxation rate V can be calculated from (22) by using 
dressed spectral densities instead of (20). This is always justified when diagrams as shown 

16 The hi<pXi/4:\ interaction in (2) also allows for scatterings without intermediate state, which we will 
discuss in the following section 6. 

17 This terminology is somewhat misleading; what breaks down here is only the loop expansion, which is 
not identical to the perturbative expansion and not controlled by a small parameter. Perturbation theory 
still works if all quantities are consistently calculated at a given order in a small parameter, which may 
require resummation [52]. The fact that such calculations are often referred to as "non-perturbative" is 
similarly misleading, as they are still based on a (resummed) perturbative expansion in a small parameter. 
This is e.g. in contrast to lattice computations in strongly coupled systems, where no small expansion 
parameter exists. 



16 




Figure 6: Diagrammatic representation of (30), i.e. (22) evaluated with the full (resummed) 
spectral density; solid lines represent (f>, dashed lines xi> dotted lines \2, the gray blobs 
symbolise self energy insertions. 



in figure 6 dominate over vertex corrections and "ladder diagrams" that contribute to 
H~(oj), which are not included in (22). For simplicity, we have constructed the Lagrangian 
(21) to suppress corrections to the g vertex and ladder diagrams. In this case the effect of 
scattering processes encoded in the diagrams shown in figure 7 can be parametrised by the 
thermal widths r, of he Xi- 

We adopt the Breit-Wigner approximation (19) for one of the spectral densities in (15) 
to perform the po integration. For the other spectral density we need to keep the full u 
dependence of Il~ in the numerator to account for pfp nt (po), especially near the regions 
where the arguments of the Bose-Einstein distributions vanish. This corresponds to the 
approximation (18). Then 11" (w) for q = reads 



rioH 



g_ 

2 



d\y\ 

(2ny 



(Z 2 U^(oo - fix) - e2t(co - fii)) (1 + / B (n0 + f B (u - fix)) 



->Z 2 W 2 (u — fix) + e(u — Q\ 



(Z 2 U 2 (u + fix) - e2i{u + fix)) (f B (u + fix) - M^i)) 



SI! (w + 



*\2 



] -,z 2 n 2 -(u + n* 1 ) + e(Lj + n 1 ] 



(Z 1 U^(co + n 2 ) - e2i{u + 2 )) (/ B (w + o 2 ) - f B (n 2 )) 



'cu + n 2 y 



iz 1 n~(u} + Ci 2 ) + e(uj + n 2 ] 



(z^icu - n 2 ) - e2i(uj - n 2 )) (i + f B (n 2 ) + f B (u - n 2 )) 



n* 2 



(30) 



Here we have again suppressed the spacial momentum index, f2 f = Q ip and Tj = T ip . The 
Xi self energies Il~ are all evaluated for spacial momenta p. The infinitesimal parameter 
e has to be kept to regularize the integrand in case a pole lies outside the support of the 
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Figure 7: Finite Ti contributions to (30) for Xi interactions as in (21) in diagram a) and for 
Xi interactions as in (32) in diagram b) . Solid lines represent <f>, dashed lines xi > dotted lines 
X2, double lines £. 

in the numerator 18 . We kept the factors Zi, which are in good approximation unity 
for the quartic self-interaction in (21), because we will need them in the following section. 
The result holds to first order in Ti/Qj. In the limit Ti/Qi — > it reproduces the analytic 
formula (25). In contrast to (25), it does not vanish for T c < T < T c because it includes 
processes in which the Xi are off-shell. Mathematically it is easy to see that the RHS of (30) 
does not vanish: The integrand is positive semidefinite, certainly non-vanishing in some 
regions and integrated over the entire positive |p|-axis. Physically this means that there 
exists no strictly stable (quasi)particles at finite temperature because even those particles 
that would be stable in vacuum can be annihilated in a scattering with quanta from the 
bath. 19 

4.3 the kinematic regimes 

We consider the case A 2 ^> Ai to study the temperature regime T > T c . This e.g. qualita- 
tively resembles a situation in which one of the fields is charged under some gauge group 
while the other one is not. In the SM a similar situation is realised for the Higgs coupling 
to electrons and neutrinos. It also resembles a coupling to active and sterile neutrinos. 
Then off-shell transport due to the tiny ^-width is negligible (while we keep the % 2 width) 
and we can use the functions T>q\ Sq^ introduced in (25) to approximate (30) by 

u-(M) 9 - f M?! ( n 2 -(M - QQ (1 + /gffji) + f B (M - go) 

° l } 2j(2nyn 1 { ((M — fii) 2 — Q%) 2 + (r 2 fi 2 ) 2 

n 2 -(M + go (mm + np - MgQ) \ 
((M + no 2 -fi|) 2 + (r 2 n 2 ) 2 J 

+i{V [ £\M) + S^{M) - V [ 9] (-M) - S [ 9] (-M)), (31) 

18 The support of the full spectral self energy covers the entire energy axis at T ^ because it contains 
cuts corresponding to scattering processes without kinematic branchcut. However, at a given order in 
perturbation theory, the support may be limited due to such branchcuts, cf. (25), (33) and (65). 

19 Similar behaviour has been found in the thermal production rate of right handed neutrinos at high 
temperature [50], which is relevant for a consistent treatment of leptogenesis, and chemical equilibration 
in hot plasmas [48]. 
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The integral in (31) has to be evaluated numerically. This requires knowledge of the self 
energy 11^ for arbitrary momenta and off-shell energies. We obtain these by creating a 
numerical table, based on the equations given in appendix A, which we then insert into 
(31). The result for T obtained from (31) is shown and compared to the on-shell result 
from (25) in figure 8. 

The behaviour of (25) can be understood easily. For T < T c dissipates its energy via 
decays and inverse decays. The temperature T max at which this is most efficient is deter- 
mined by the competition between amplification by Bose enhancement and the decreasing 
phase space at increasing temperature. For T > T c the decay is kinematically forbidden and 
(25) is zero. For negligible vacuum masses, T c can be estimated as T c ~ my/24:/ (y/Xi+y/X^). 
At T c ~ m\/2A/\\f\i — \f\2\ the processes 0Xi X2 become kinematically allowed and 
give rise to Landau damping. 

The full relaxation rate, calculated from (31), shows a very similar behaviour. It also 
increases linearly with T due to the effect of Bose enhancement between T and a tem- 
perature T' max ~ T max , where (31) has a local maximum. Beyond T' max the relaxation 
rate again decreases with T, as the shrinking phase space overpowers the effect of Bose 
enhancement. In the region T c < T < T c , where at least one of the intermediate states in 
any process contributing to the relaxation, we observe T 7^ due to off-shell transport. T 
receives contributions from scatterings encoded in the loop diagram shown in figure 3. As 
an example we show one scattering process (0Xi ^ X2X2X2) obtained from the cut in figure 
3, but the loop also contains scatterings 0X2 X1X2X2, 0X2X2 X1X2 and 0XiX2 X2X2 
as well as (if kinematically allowed) decays and inverse decays such as o X1X2X2X2 
or xi "H- 0x2X2X2- These are suppressed by an additional coupling A2, the virtuality of 
the intermediate X2 quasiparticle and the required phase space. They are not included in 
(25). On the other hand, they are amplified by the infrared behaviour of the Jb due to 
Bose enhancement. The resulting overall suppression is very efficient not only compared 
to T(T max ), but even relative to the decay width in vacuum. Thus, if this scenario is real- 
ized in nature, perturbative reheating indeed effectively stopped at T c , as claimed in [34] 
(the plasma would, however, be heated much faster than calculated in that article as the 
authors neglected the enhancement by Bose enhancement for T < T c ). 

The effective suppression is a consequence of the interaction Lagrangian (21). The 
suppression in the regime T c < T < T c is particularly strong for the interaction Lagrangian 
(21) for three reasons. First, it does not contain any interaction that couples Xi an d X2 
to each other. 20 If xi an d X2 were e.g. charged under the same gauge group, T would 
receive contributions from "ladder diagrams" similar to the one for fermions shown in 
figure 16. We avoided such terms to obtain the transparent expression (30), which is 
analytic up to a one-dimensional integral. It is known from other examples that such 
contributions can qualitatively change the behaviour. For instance, the production rate 
of right handed neutrinos from a thermal plasma at tree level exhibits a kinematically 
forbidden temperature region similar to the interval T c < T < T c discussed here [49]. 

20 We indirectly assumed that they communicate in some way via £bath by assigning the same temper- 
ature T to both. The results in this section remain valid as long as the kinetic equilibration between xi 
and X2 happens faster than 1/T. 
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There ladder diagrams do contribute, and in [50] it has been found that their resummation 
can almost entirely overcome the suppression in the "forbidden region" . For a choice of Cx 
that yields contributions to V from ladder diagrams during reheating we expect a similar 
effect, though it is probably less pronounced than in the neutrino production rate due to the 
absence of the collinear kinematics that is crucial to the result found in [50] (cf. caption 
of figure 16). Second, the quartic self interaction in (21) leads to a contribution to the 
dispersive part of the Xrself energy at order Xi, cf. figure 5a), giving rise to a thermal mass 
oc y/Xj,T. In contrast to that, the lowest order contribution to the dissipative self energy Il~ 
is of order A?, cf. figure 5b). This disfavours off-shell transport if Aj is a small parameter: 
Il~ appears in the numerator of (30) and scales as oc Af. The critical temperature T c , on 
the other hand, depends on the coupling as A^ 2 because Mi oc X]^ 2 . The factors fs in the 
numerator of (30) enhance T in the "forbidden region" more efficiently for a large T c , thus 
small Aj. This effect competes with the X 2 suppression from Ilr in the numerator. For a 
gauge or Yukawa coupling the dispersive and dissipative part of the self energy appear at 
the same order, and the thermal mass is proportional to the coupling (rather than its square 
root). This allows the effect of Bose enhancement to be more efficient at T c . Furthermore, 
the fact that the lowest order contribution to Il~ comes from a two loop diagram leads to 
an additional suppression, which e.g. can be noticed in the small numerical prefactor in 
(40). This furthermore suppresses T in the forbidden region via the II - in the denominator 
of (30) 21 . Finally, the momentum independence of the thermal mass (24) implies that 
even soft modes become very thermally heavy. For these reasons the upper bound on the 
reheating temperature found in this section is a rather special case. 



5 Example II: no blocking by thermal masses 

We now consider a situation in which the coupling between and Xi is the same as in 
(21), but the Xi have different interactions. This illustrates how strongly the behaviour of 
T depends on the interactions within the primordial plasma even if the inflaton coupling 
is the same. We use the following model 

£ = \d»<pd^- l -m 2 4> 2 - g<p X iX2 (32) 
+ Yl (l^Xid^Xi - \rnlx] + - 0i rf) + ^bath- 

i=l ^ ' 

Here £ is another real scalar field. We can again use (22) to calculate T, but in this case 
the dissipative and dispersive parts of the x% se ^ energies are both calculated from the one 

21 A small IT^ does not suppress T in the regions where on-shcll transport is kinematically allowed 
because in this case the integration volume includes narrow pole regions where the denominator of the 
integrand in (30) is of order 1 / ) 2 and the integral is of order one. 
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Figure 8: The rate T in the model defined by (21 ) as a function of T, normalized to its 
zero temperature value. We set Ai = 0.01, A2 = 1. mi = m-2 = 0.001m. The solid line is 
calculated from (25), it takes into account thermal mass corrections for Xi, but neglects the 
thermal widths. The dashed light blue is the result one obtains when also neglecting thermal 
masses. The dark blue dots are numerical evaluations of (31 ), which also take into account 
thermal \i-widths. 



Figure 9: Leading order contributions to the Xi- se ^~ enei S^ es m the model (32); we use the 
same notation as in hgure 7. 
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loop diagrams shown in figure 9. The dissipative part can be found analytically and reads 



Ib{uj ) f B (u + -Pq) 



9{p 2 -{2m i ) 2 ) [p ^-p ] (33) 



with 



2 2 y p 2 K J 

Since the dissipative and dispersive self energies can be identified with real and imaginary 
part of the function ITf , one can use (11) and the Kramers- Kronig relations to find the 
following equation for Q 22 

(35) 

The dispersive part is plotted in figure 10. The most prominent feature of the function 
Rellp(po) is the spike, which appears around p\ = p 2 + (2m^) 2 because the pole in the 
last term in (35) passes the threshold of (33). The momentum dependence implies that 
the solution to (35) in this regime can in general not be parametrized by a momentum 
independent thermal mass. It is sometimes argued that one can simply use an effective 
thermal mass that is determined by fixing the momentum to T, as most particles in a 
thermal plasma at T > have momenta ~ T. This argument fails in our case because 
the integral in (30) is not always dominated by the region around |p| ~ T that maximizes 
P 2 /b(|p|), but by the regions around the poles. The momentum scale in <fi — > XiXj decays 
is set by the <p mass, the maximal momentum the decay products can have is |p| = m/2. 
For T ^> m this corresponds to "soft momenta". That means that, for sufficiently high 
T, T is sensitive to the shape of the Xi spectral densities in the infrared regime, where 
the principal value in (35) dominates over p 2 + mf. In this regime the pole structure can 
be rather complicated. Physically this happens because particles with soft momenta can 
strongly affected by the medium. 

The non-monotonic behaviour shown in figure 10 makes it possible that (35) has more 
than one positive solution. The additional solutions could be interpreted as a collective 
excitation. On the other hand, there is no guarantee that (35) has a solution at all because 
the principal value tends to be negative for small p. This suggests that for some parameter 
choices there exists a minimum momentum for quasiparticles. For smaller momenta the 
spectral density (9) is smooth, with no poles or narrow Breit-Wigner peaks, and there is 
no propagating excitation with a quasiparticle interpretation. Physically this would mean 
that the thermal damping in the plasma prevents excitations with small momenta from 
propagating; the resonances are so broad that their lifetime is comparable or shorter than 



22 When evaluating (35) we drop the T = part of (33), which would give rise to an infinite contribution 
that has to be removed by mass renormalization. As for (24), we assume that this term has already been 
absorbed into m*. 
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their frequency . Whether or not a well-defined resonance can be identified depends on 
the value of the RHS at its minimum. 24 We restrict the analysis to the case where (35) has 
exactly one solution, for which the hierarchy Tj <C fi, holds and allows to define a (possibly 
finite width) quasiparticle. Then (30) can be used to calculate T. As a cross-check for the 
validity of (30), we evaluated the sum rule 



for each point in the integration volume for each temperature. It directly follows from 
the commutation relations \x(xi), x( x 2)]\t 1 =t 2 — iS(xi — x 2 ) for a scalar field. Using the 
sum rule we find that the approximations (19) and/or (18) can fail at high temperature 
(T rrii,m) for very hard (|p| ^> m;) and very soft (|p| <C m^) momenta. This can easily 
be understood qualitatively. For soft momenta the principal value dominates (35). As 
discussed before, in this regime the pole structure can be complicated and (18), (19) may 
not provide useful approximations. For hard momenta one can in good approximation 
neglect all vacuum masses. The thermal correction is smaller than jjjT, hence the pole 
of pi and the spike- feature in ReII~ (cf. figure 10) approach each other near p$ ~ |p|. 
The underlying assumption of the Breit-Wigner approximation is that U R does not change 
considerably across the quasiparticle peak and can be replaced by its value at Qi in loop 
integrals. But near the spike, Rellf as a function of po is steep and the factor Zi is 
considerably different from zero. Physically this means that the pcrintegral in (36) is not 
strongly dominated by the peak region and the continuum part of pi cannot be neglected. 
In this case (18) may still hold, but (19) fails . In the numerical results plotted in figure 
11 we chose parameters such that the deviation from (36) is less than 15% in the regions 
that contribute significantly to the integral in (30). 

The thermal mass shift depends on both, temperature and momentum. It can have 
either sign, dependent on the choice of m,, T, Qi and p. It is therefore impossible to 
define a unique temperature T c at which the decay <fi — > X1X2 becomes kinematically 
forbidden. It can be forbidden for some Xr m odes, which have received a sufficiently large 
positive "thermal mass" at a given temperature, but is still allowed for other modes if these 
have received a smaller (possibly even negative) correction to the dispersion relation. We 
therefore always have to insert the full dispersion relation into (30) 25 . Figure 11 shows 

23 0f course, the non-existence of a solution for (35) for a given mode p does not necessarily imply that 
there can be no quasiparticle with that momentum: (9) also exhibits a pronounced peak if the LHS is 
sufficiently small for a po > 2mf , i.e. above the Xi ~^ threshold. There appears to be no quasiparticle 
if the LHS of (35) is small, but non-zero, for \po\ < 2m^ because in this region only exact zeros of the 
denominator give a finite contribution to the spectral density (9) in the limit e — > 0. However, this 
statement has to be taken with care because the resummed H~ in the denominator of pi is non-zero and 
may be unsuppressed in the regime where the decay x% ~~ ^ is forbidden, in the same way as we illustrate 
it in this section for IT - . In our numerical analysis we used the one- loop result (33) for simplicity. 

24 Note that in this case (35) and (15) cannot be used to compute the width and precise dispersion 
relation, but the denominator of (9) has to be evaluated in detail. 

25 Strictly speaking we cannot even apply (30) because in its derivation we assumed = but it 
should still hold as long as M\ <C ff. 
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Figure 10: The dispersive part of the self energy IT^po) = Hf (po)/<?f at T = m for \p\ = 
0.1T (solid line) and \p\ = T (dashed line) as a function of energy, with negligible ^-masses 
m s /T = 10~ 5 . 



the temperature dependence of T as a function of T. For simplicity, we neglected thermal 
corrections to m^. We checked that for the parameters we plot, a generic thermal mass 
does not change the results qualitatively. There is clearly no suppression due to thermal 
masses even for T 3> m,mj. This is relatively generic if one chooses rrii and Qi less than an 
order of magnitude smaller than m. The reason is that the decay products have momenta 
< m/2. For rrii ^ m <^ T, these are soft, i.e. much smaller than the typical momentum 
~ T in the plasma. For soft momenta the correction to the Xi dispersion relations are small 
or even negative. That means the dispersion relations for very soft Xr m odes are almost 
free of thermal corrections and for these modes there are no large thermal masses even at 
high temperature. 26 The high sensitivity of the reheating rate to the dispersion relations 
at soft momenta is due to the fact that the inflaton decays at rest, hence the momenta of 
the decay products are determined by m and not of the scale T. This is in contrast to most 
other transport phenomena in thermal equilibrium, where most of the energy is stored in 
modes with hard momenta ~ T which receive thermal mass corrections of order aT. 

26 The possibility of a vanishing thermal mass has been studied in a different context in [72], where 
indications were found that the thermal mass may vanish in strongly coupled gauge theories. 
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Figure 11: The rate T in the model (32) as a function ofT, normalized to its zero temperature 
value. We set 0i = 02 = 0.1m, mi = mi = 0.2m. The dashed light blue is the result one 
obtains when also neglecting thermal masses. The dark blue dots arc numerical evaluations 
of (30), using the dispersion relations obtained from (35). 



6 Example III: reheating by Landau damping 

We now consider the case that the Xi se lf interactions are the same as in (21), but their 
coupling to is different. We use the following model 

i=i ^ • • / 

(37) 

The leading order contributions from the quartic interactions to II - come from diagrams 
of the type shown in figure 2b) and c). In analogy to (22), they can be calculated from 

x [(1 + f B (po)) (1 + f B (ko)) (1 + f B (!o)) ~ fB(Po)fB(k )f B (l )] . (38) 
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The integrand is again given by a product of spectral densities, weighted with distribution 
functions. In the on-shell limit (20) it reads 



X ( (1 + / p ) (1 + /k) (1 + /l) - /p/k/l) 

+3(/p(l + / k )(l + /i)-(l + / P )/ k /i)) 

+ fi ip - fi ik - Qii) - 8{u - Q ip + fi ik + fia)) , (39) 

with f p = fsi^ip) etc. Again the integral is only non- vanishing when the submanifolds 
defined by the on-shell (^-functions intersect somewhere in the integration volume. In 
vacuum, this would only be possible as long as the decay (f> — > XiXiXi is kinematically 
allowed, i.e. for M > 3Mj. However, at finite temperature (38) also includes terms that 
can be interpreted as scatterings (pXi ^ XiXi- They give a contribution to T that allows 
to reheat the plasma when the decay <p — > XiXiXi is forbidden. In contrast to the trilinear 
interaction discussed in the previous section, these processes do not require an intermediate 
particle. This is a generic feature of vertices that connect more than three propagators. 

However, since scatterings require Xi quanta in the initial state, they are suppressed at 
low temperatures, when there are not sufficiently many scattering partners. In the high 
temperature regime (39) can be calculated analytically. For our purpose, the most relevant 
case is T ~ T c . In [65-67] an analytic expression has been obtained for M = Mj, 2 



E 



h 2 T 2 

for M = Mi. (40) 
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This result should be approximately correct in the regime T ~ T c , where u = M = 2Mj. 
It is also interesting to consider the case Mi 3> M. In this regime we find, to leading order 
in M/Mi, 

r^JpLll( 1 + l0 JZAMi)) t0lM «M,. (4i) 

6(2vr) 4 M?\ B \8MJJ K J 

In figure 12 we plot T as a function of T from numerical evaluation of (39). The analytic 
results (40) and (41) provide rather accurate approximations in some temperature regimes. 
(39) only takes into account tree level processes (ft <H- XiXiXi ari d (pXi XiXi with all 
involved quasiparticles on-shell. In order to incorporate processes with more particles in 
the initial and final state as well as intermediate off-shell quasiparticles, one would have to 
include the finite widths, which are neglected in (20), and use resummed vertices. 

Let us combine the models (21) and (37). We assume that and Xi couple by the 
trilinear g<pXiX2 as well as the quartic /ii/(4!)0xf interactions. We compare the efficiency 
of decays -H- x 1X2 to Landau damping by scatterings 4>Xi -H- XiXi i n the regime T < T c , 
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Note that the definition of T varies by a factor 2 in different articles. 
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Figure 12: The contribution to T from a /ij/(4!)0%? -interaction for Xi = 0.1. The solid line 
is a numerical evaluation of (39), the dashed line the approximation (40) and the dotted line 
the approximation (41). All contributions are normalized to the value at T = 0. For larger 
temperatures than plotted here our numerical evaluation is not reliable, but T may decrease 
again due to the shrinking phase space for <pXi -H- XiXi scatterings. Note also that (39) does 
not take into account hnite widths or resummed vertices. 

both on-shell. 28 We consider the simplest case, with mi = m 2 = m x , h\ = hi = h and 
Ai = A 2 = A. The rate (25) then reduces to 

1/2 

{l + 2f B (M/2))9(u;-2M x ), (42) 

The rate for cf>Xi XiXi scatterings can be estimated by (40). At T c « (6(m 2 — 4m 2 )/A) 1 / 2 
it takes the value 

, mS h 2 m ( /2m v \ 2 \ . , 

r^(T c )«2— (l-(-^)j. (43) 

This can be compared to the decay rate (42) at the temperature where it is maximal and 
at T = 0. The temperature dependence of (42) is plotted in figure 13, below T c it is very 
similar to the Mi ^ M 2 case discussed in section 4.3. For < T < T c <p dissipates its energy 
via the decay <p X\Xi- Since both final states are bosonic, Bose enhancement causes V 
to increase as a function of T. At the same time, the phase space for this process shrinks 
due to the large thermal Xi masses. Competition between these two effects determines 
the temperature T max , at which T is maximal. The value of T max can be found using the 
requirement dY& ecay /dT = at T = T max , which allows to formulate the condition 

m (6((2m x ) 2 - m 2 ) + T 2 ax A) + T^ ax Xsmh(m/T max ) = 0. (44) 

28 The discussion closely follows [38]. 
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Figure 13: The decay rate (42) for mi = m,2 = 0.1m and Ai = A2 = 0.1. 



For T max ^> m the sinh can be Taylor expanded to obtain 



I (m 2 - (2m x ) 2 ) 



1/2 



With M 2 (T max ) « (m 2 + (m/2) 2 )/2 this allows to estimate 



decay (^mai 



<:r 



2m x N 
m 



87rm V A 

where T maa; 3> m allowed to expand /s(m) in m/T. Comparison yields 



Tdecay (T max ) 
F decay (^ = 0) 

Tscatter (-^c) 
F decay (^ = 0) 

r scatter (-^c) 
Tdecay y^max ) 




A 
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2m x x 
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1/2 



2m x x 
m 



1/2 




scatter V- 1 c 



24A V 24r decay (T = 0) 



m v =0 



(45) 



(46) 



(47) 
(48) 
(49) 



Remarkably (49) does not depend on m x . In the interesting case m x <C m also the other 
ratios are very simple. The parameter that governs the maximal amount of amplification 
of the decay rate by Bose enhancement is ~ (24/A) 1 / 2 . The parameter that determines 
whether relaxation via scatterings is efficient at T c is (h 2 /X) ■ (m/g) 2 . For g/m ~ h the 
rate r scatter (T c ) can easily be bigger than r decay (T = 0) or even T dccay (T max ). 
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7 Example IV: production of fermions 



Matter in the real world is composed of fermions with gauge interactions. So far we have 
only considered the case that first dumps its energy into other bosons and assumed that 
fermions get created in the subsequent decay chain or inelastic scatterings in the plasma. 
This is justified because the transition into bosonic final states is usually Bose-enhanced, 
while those into fermionic final states are Pauli suppressed. In this last example we study 
T in a model where directly couples to fermions, 

£ = l -d^(j> - l -m 2 <f + * (i? - m) - Y<pm - oM-fAJS! - + £ bath . 

(50) 

Here is a gauge field and is the field strength tensor. The rate V calculated from 
the diagram shown in figure 2d) is given by 

— ^tr(p p (p )Pp- q (po-u;)) (51) 
x ((1 - /f(Po)) (1 - /f(w - Po)) - /f(Po)/f(w - Po)) , 
where /f(u;) = [& } I T + l)" 1 is the Fermi-Dirac distribution. 



7.1 Temperatures T <T C 

For T < m thermal corrections to the fermion properties are negligible and we can use the 
bare fermion spectral density to evaluate (51) 

p[T(R)) = 27rsign(p )(y+ m)5(p 2 - m 2 ). (52) 

The result is 

3/2 

[1 - 2f F (m/2)} 0{m - 2m) for T < m (53) 

In the more interesting regime T 3> m (we neglect m in the following) we again have to 
use resummed spectral densities as in the scalar case. In the hard thermal loop (HTL) 
approximation one obtains [69] 

Pp(Po) = 2 ((7o - P7)P+ + (7o + P7)p-) • (54) 
Here p>7 = Pi7t/|p|. The functions 

p ± (p) ~ 2n[pl olc (p) + p c ± ont (p)] (55) 



r Y2 
1 = — m 

8tt 
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are the sum of singular contributions p± ole and continua p™ nt . The poles of the singular 
parts define quasiparticles with energies Q.±, 

Pl° lc (p) = Z ± 6( Po - £1 ± ] 
the continuous part is given by 



Z T 6(p + Q. 



(56) 



p^(p)=0(l-x 2 )f-(lT 
P 



x I 



X 



1 =F x ± y 2 ( (1 =F x) In 



X + 1 




x — 1 


«) 





TrY (1 =F x) 2 



(57) 



We used the notations x = po/|p| and y 2 = 9Jt 2 /(2p) 2 . The thermal mass reads 971/ = 
cxT/2. 9Jlf is sometimes referred to as asymptotic mass, it differs by a factor \/2 from the 
plasma frequency CO/ at |p| = 0. The residues are 



4y2p2 

The dispersion relations Cl + and — Q_ are the solutions to 

x + 1 







IPI 



1 + y 2 (1 - x) In 



(5f 



(59) 



There are two types of quasiparticles, with dispersion relations Q + and D. . The former 
can be interpreted as one-particle states that are screened by the plasma, the latter are 
collective resonances often dubbed "holes" or "plasminos". This complicated structure 
makes it obvious that correct results for fermions cannot be obtained when simply replacing 
intrinsic by thermal masses. The dispersion relations D.± can be expressed in terms of the 
Lambert W-function as [55] 



a. 



IpI 



W-i(s) 



a 



W-x{s) + 1 
-u/+i)_ They are plotted in figure 14. 



iPl 



Wn(s) 



W (s) + 1 



(60) 



with s = e 

At temperatures T <C M/cc we can approximate (51) by neglecting p^? nt . Since m ^> 
m, the daughter particles have hard momenta |p| ^> m, OJt/. In this regime 29 one can 
approximate 



y 2 (l + 21n(y)) 



a 1 y 2 exp 



-y~ 



and 



0+ ~ |p|(l + 2y 2 ) ~ (p 2 + 9JT 2 f ) 1/2 Q_ ~ |p|(l + 2a- 1 exp(-y- 2 )) 



(61) 



(62) 



29 In principal the HTL approximation only holds for soft momenta |p| -C T. We here use it for arbitrarily 
hard momenta, assuming that medium related modifications of the spectral density play no role for |p| > T. 
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Figure 14: The fermionic dispersion relations Q+ (red line) and O- (blue line) in units of the 
plasma frequency Wf = DJlf/\/2. The dashed line shows the approximation (p 2 + DJl 2 ?) 1 ^ 2 for 
hard momenta. For simplicity we assumed in (50) that have electromagnetic interactions; 
the behaviour in a QCD plasma is similar. 

The dispersion relation Q + for the screened particle approaches that of a particle with 
mass ~ DJtf = otT/2. This thermal mass only appears in the dispersion relation and does 
not break the chiral symmetry [69], as one can see from (54). The holes effectively become 
massless for hard momenta ^> 97t/. However, the residue Z_ is exponentially suppressed 
for large momenta. The physical reason is that holes are collective phenomena. Such 
excitations appear at length scales which are associated with soft momentum modes ~ <xT 
and longer than the typical inter-particle distance ~ 1/T. The occupation numbers of soft 
modes ~ <xT in a thermal plasma are small compared to the hard modes ~ T. Therefore 
it is usually argued that processes involving holes do not contribute significantly to the 
damping unless other channels are forbidden or they are resonantly amplified 30 . Thus, the 
spectral density can be approximated as 

Pp(Po) - 7T (( 7o - P7)<Kpo - + (7o + f>j)5(p + 0+)) . (63) 
Insertion of (63) into (51) yields for |q| = 

1/2 

(1 - 2f F (M/2)) 6(M - 2M f ) for m < Wl f < M (64) 

where we approximate + ~ (p 2 + VJtj) 1 ^ 2 . The analytic formula (64) for small a is in 
good agreement with the numerical results found in [54] and our own results shown in 

30 Similar arguments apply to hydrodynamic and ultrasoft modes studied in [70]. 
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Figure 15: Cuts through the <p self energy with resummed fermion propagators. 

figure 17. Note that it cannot be obtained from (64) by simply replacing the bare mass 
m with the thermal mass 971/. We can draw two important conclusions in the regime 
9JT/ < m < T. First, if the 0-couplings to fermions and bosons are of comparable strength, 
then comparison with (42) shows that the contributions to V from decay into fermions 
is suppressed due to Pauli blocking. Second, the effect of the thermal mass 97t/ on T is 
relevant only for T > m/cc, while Pauli suppression is effective already at T > m. 

7.2 Temperatures T > T c 

It is, however, not clear that the above conclusions hold at higher temperature. The reason 
for the Pauli suppression for T > m is that the energy of the \l/-quasiparticles in <fi — > \]/\]/ 
decays is ~ m/2, independently of the temperature. Cuts through the resummed prop- 
agators as shown in figure 15 include scatterings -H- Wj. In these Landau-damping 
processes, the fermions can have hard momenta and energies ~ T, hence there is no Pauli 
blocking. The continuous parts p™ nt of the fermion spectral densities consistently incorpo- 
rate the leading order contribution from Landau damping in the resummed propagators. 
They do, however, not include contributions from "ladder diagrams" as shown in figure 16. 
Inclusion of these would require not only resummed propagators, but also resummed ver- 
tices. In [73] it was argued that corrections to the decay rate of a scalar into fermions with 
gauge interactions from resummation of vertices are of higher order, see also [74] 31 . This 
implies that evaluation of (51) with the full resummed spectral densities (54) is sufficient 
to include all processes at leading order in oc. We split the product of spectral densities in 
(51) in three parts. 

The first part is related to the red cut in figure 15 and contains only contributions from 

31 Note that there is no Ward-Takahashi identity for the Yukawa vertex [74], as there would be for gauge 
interactions. 
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Figure 16: A ladder diagram with one rung. For the Yukawa coupling considered here ver- 
tex corrections from such diagrams are sub-dominant [74]. It has been found in [50] that 
rcsummation of ladder diagrams can give a large contribution to relaxation rates, cf. also 
[75] in that context. However, in that work the momenta of typical external particles is ~ T. 
Therefore the particles have collinear momenta in the rest frame of the bath and propagate 
parallel for a sufficiently long time to exchange several gauge bosons. It is precisely this 
collinear enhancement that leads to the large contribution. Furthermore, one of the internal 
lines in that calculation is a boson. We therefore expect that the effect of ladder diagrams is 
smaller in the case of reheating and does not change our hndings qualitatively. 

the on-shell parts p± ole of the spectral densities. It reads 



n (™)lpx P = -^ J d\p\ P 2 



Z\ [1 - f F (Cl + ) - f F (m - Cl+)]5(m - 20+) 
+ Z + Z_ [1 - f F (Cl + ) - f F (m -CL + + - f F (Cl_ + m)]5(m + D_ - D_ 

+ Z 2 _[l-f F (Cl^-f F (m-CL^]6(m-2Cl^ 



(65) 



Since the finite widths are neglected in (56) just as in (20) for the scalar case, evaluation of 
(51) with (56) leads to a qualitatively similar behaviour as in (25); in particular, tempera- 
tures T c and T c can be defined. The ^-function allows to evaluate the integral analytically, 
but the zeros of their arguments have to be found numerically. The contribution to V from 
(65) is shown as red line in figure 17. The term (65) contains three types of contribu- 
tions. The terms oc S(m — 2D + ) describe 0-decays into dressed \P particles. They clearly 
dominate at T < T c ~ m/oc, where (64) provides an excellent approximation. Terms 
oc 5(m — 2D_) describe the decay into two collective fermionic excitations. Since Z_ is 
suppressed for hard momenta, they only contribute significantly near T c , where the decay 
products' momenta are soft. Their main effect, visible in figure 17, is that (65) is non-zero 
for temperatures slightly above T c , where the decay into dressed particles is already kine- 
matically forbidden (2Q + > m), but the decay into collective excitations is still allowed 
(2£1_ < m). However, they also become kinematically forbidden at only slightly higher 
temperatures. Both of these contributions are suppressed by Pauli blocking at T > m. 
The contribution oc 5(m + Q„ — Q + ) comes from processes in which a <p particle and a col- 
lective \l/ excitation form a ^ particle. These become kinematically allowed above a second 
critical temperature T c , where D. + > m + Q_. 32 The contribution from such processes is 
generally small, only near the threshold at T c it is amplified by a van Hove singularity. This 



:S2 



While for scalars it was necessary that there are two fields Xi with Mi ^ M 2 for the similar processes 
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singularity originates from the infinite density of states (dp)/(dD_) at the local minimum 
of 0_ for |p| ^ visible in figure 14. For most parameter choices physical observables are 
not affected by the singularity because the density of states only diverges in one point and 
gives a finite contribution when convolved with smooth functions. Only if the mass shell 
defined by the zero of 5(m + Q_ — Q + ) happens to coincide with the van Hove singularity 
T is divergent. This happens at T c . The divergence is, however, not physical and can be 
removed by taking into account the finite width of the quasiparticles, which is neglected 
in the hard thermal loop approximation (56). This simply means that we cannot trust 
our approximations in a small temperature interval around T c , which we exclude from the 
following considerations. Similar behaviour was found for the dilepton production rate in 
a quark-gluon plasma [79]. 

The second term comes from cuts through the 0-self energy that cut through one \l/-self 
energy (blue in figure 15) and reads 



n M|px C 



-iY 2 ^ j rf|p|p 2 



(i - f F (a + ) - f F (m - n + ))p cont (n + - m ) 



+(Mo+) - Mn + + m)) P c _ ont (o + + 



m) 



(1 - MO_) - f F (m - Q_)) p c + ont (n_ 



+ (/ F (n_) - MO. + m)) p c + ont (0_ + m) 



(66) 



The continuum parts p± nt originate from the imaginary parts of the \l/-self energies, hence 
including them means that we take into account dissipative effects for \1/ though the thermal 
width is neglected in (56). It includes processes such as -H- 7$, with intermediate 
These are kinematically always allowed, hence (66) is nonzero for arbitrary temperatures 
there is no equivalent to T c . However, they always contain fermions in the initial and final 
state. Since the momentum of the outgoing fermion kinematically must be smaller than 
that of the incoming one (because the is at rest), they tend to be Pauli suppressed. 
Indeed this contribution never exceeds the vacuum T in figure 17. 
The last term reads 



n (m)\cxc 



x 



-zF 2 ^ J d\p\dp p 2 [l - f F (po) - /f(t» - Po 
[p c + ont (p )p c _ ont (Po - m) + p c „ ont (p )p c + ont (Po - m)] 



(67) 



It contains the contributions from cuts through the 0-self energy that cut through both 
self energies (green cut in figure 15). They can have more particles in the final state. 
This increased flexibility in the phase space leaves more room for the fermions to have 



fiXi Xj to be kinematically allowed, they are always possible for fermions because the hole provides a 
second particle with different "mass" . 
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hard momenta and avoid Pauli suppression. At the same time, transitions with bosonic 
final states with soft momenta are enhanced. At low temperatures such processes are 
negligible because they are of higher order in the coupling, but at high temperature this 
suppression is cancelled by the high density of scattering partners (formally the Bose- 
Einstein distribution). At temperatures T 3> T c the term (67) becomes the dominant 
source of dissipation; it continues to grow with T and can exceed the vacuum decay rate 
by orders of magnitude. 

While it is usually assumed that fermionic channels are not important for reheating 
because they are Pauli suppressed, our result suggests that they can give a huge (possibly 
dominant) contribution to T at very large temperatures. If the Yukawa coupling is the only 
interaction that has, then the universe will never reach this temperature regime because 
of the suppression of T around T < m/oc. But if some other interaction is efficient in this 
regime, then processes involving fermions can be crucial at higher temperatures. 

8 Conclusions 

We calculated the relaxation rate T of a scalar field in a plasma of scalars with trilinear 
and quartic couplings and fermions with gauge interaction using resummed perturbation 
theory. We focussed on the regime where a quasiparticle description applies and the dissi- 
pation is dominated by perturbative processes involving individual quanta rather than par- 
ticle production from the time varying background field (0) . We paid particular attention 
to the temperature dependence of the phase space, caused by the modified quasiparticle 
dispersion relations in the plasma, and to quantum statistical effects. 

Many transport phenomena can in good approximation be understood in the descrip- 
tion of a hot plasma as a collection of quasiparticles with momentum independent thermal 
masses. The reason is that for weakly coupled theories at temperature T in an adiabati- 
cally changing background, the dispersion relations for quasiparticles with hard momenta 
p ~ T asymptotically approach u ~ (p 2 + M 2 ) 1 / 2 , where Mj ~ aT is a momentum inde- 
pendent thermal mass. Since most particles in a thermal plasma have momenta ~ T, this 
description is appropriate for most microphysical processes. 

This is, however, not the case for the decay of heavy particles with mass m at rest or 
the dissipation of the oscillating zero mode of a field, such as the inflaton. The reason is 
that in such a decay the typical momentum of the decay products is not ~ T, but of order 
m/2 or smaller. For T ^> m such momenta are soft compared to typical energies ~ T in 
the plasma. This makes T sensitive to the dispersion relations for soft modes. They can 
be complicated functions of momentum that differ for different interactions. 

r is of crucial importance to understand the thermal history of the universe after infla- 
tion, as it yields the rate of perturbative reheating. Our results allowed to investigate the 
possibility that thermal masses impose an upper bound on the temperature in the early 
universe by kinematically blocking the inflaton decay into the corresponding particles. We 
conclude that thermal masses may have a considerable effect on the efficiency of pertur- 
bative reheating. However, a strict bound on the reheating temperature from kinematic 
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Figure 17: The black line shows the total T in the model (50) calculated from (51) and 
normalized to its T = value. The colourful lines show the contributions from the individual 
terms (65)-(67). They can be identihed with the cuts in figure 15 as indicated by the colour 
coding. The red line is the contribution from term (65). In the regime T < T c ~ 10m it comes 
from the on-shcll decays of 4>-particles into two ^-particles. The non-zero contribution for 
temperatures slightly above T c originates from on-shcll decays of 4>-particles into two collective 
fcrmionic excitations. The contribution for T > T c ~ 60m originates from annihilations of a 
4>-particle by a collective '5 -excitations to form a ^-particle. The contributions from processes 
involving collective excitations are divergent near the thresholds due to a van Hove singularity 
in the density of states, see discussion in the text. In these regions the finite widths of the 
quasiparticles have to be taken into account to obtain a physically meaningful and finite 
result, which goes beyond the HTL approximation (54). We do not do that in this work and 
indicate the corresponding regions by dotted lines. The blue line is the contribution from term 
(66), which includes 4>^ — >• ^7 scatterings and is non-vanishing in the regime T c < T < T c . 
The green line is the contribution from term (67), which includes various scattering processes 
involving fcrmions and photons. They are not Pauli suppressed. At low temperature they are 
negligible because they contribute at higher order in <x. At high temperatures, however, this 
suppression is overcome by the high density of scattering partners and the Bosc-cnhanccmcnt 
for photons. 
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blocking of inflation decay can only be obtained in very special setups. In general, other 
mechanisms allow the inflaton to dissipate its energy into the primordial plasma at high 
temperature. These in particular include Landau damping by scattering. Even when large 
thermal masses require some intermediate particle to be off-shell, these can be relevant 
due to Bose or collinear enhancement. A main conclusion of our study is that it is usually 
not possible to express the reheating temperature in terms of the parameters appearing in 
the Lagrangian in a simple way even if one is only concerned with perturbative reheating. 
The reason is that reheating is not an instantaneous process, but takes place over a time 
interval ~ 1/T during which the universe expands. Hence, the temperature dependence 
of T is crucial to determine the time evolution of the effective temperature. It strongly 
depends on the interactions within the primordial plasma, and not only on the couplings 
of the inflaton. 

Our results may also be important for other phenomena in cosmology where thermal 
effects have been found to be relevant, including leptogenesis [3, 10, 49-51, 55, 75-78], the 
production of dark matter [40, 41], dissipation during warm inflation [44], axion produc- 
tion [43] or the evolution of flat directions [13, 42, 45, 60]. The effects we find can also be 
relevant in situations where perturbative and non-perturbative processes work together, 
such as combined reheating [30], curvaton decay [46] or during Affleck Dine baryogenesis 
[47]. 
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A The Xi self-energy 

In order to evaluate (30) we need IL^ (po) for off-shell energies and arbitrary momenta. 
The mass correction is dominated by the tadpole diagram shown in figure 5a) and can be 
estimated by (24). The self energy n - for Xi is given by the diagram shown in figure 5a). It 
can be obtained from (39) with the replacement hi — > A,. This integral can be conveniently 
rewritten in the form 



with 



Tl pi ( Po ) = 2i(2?W(po) + S]?(pd) - pW(-po) - S®(-p )) (K 

v p^ = % ^hs 6 bo) / dQ k / dn x F d {^n hP{h A)x{z) 

Z4{ZTV) J M . J M . 

\2 poo poo 

5 p(Po) = ^-r^—^p,) dn k dQiFsfok, Q\,p , —A)l(Z) 
24(2tt)° } Mi J n + 
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with 



F d (Q k ,n h A) = (i + /fl(ni))(i + /B(n k ))(i + /fl(A))-/B(n 1 )/fl(nk)/B(A) 
F s (n k ,n h A) = (i + f B m)(i + fB(^))fB(A)-f B (n l )f B (n k )(i + f B (A)) 

0(I(x,y,Z)) 



X{Z) = J dx J dy- 



y/I(x,y,Z) 
I(x,y,Z) = (1 - x 2 )(l - y 2 ) - (Z - xyf 



A 2 - B 2 + 2y|p||k 


+ 2x 


Pll 


1| 


2|k|| 


1| 



A = po - ft k - ft! , B z = p 2 + k' + V + Ml 

ft ± = max [po - ft k ± Mi, Mi] , Z 

The variables x nd y are the cosines of the two nontrivial integration angles. Note that 
the integration limits imply that T>p (po) vanishes for po < 3Mj, as expected from energy 
and momentum conservation in the decay and inverse decay processes %% ^ XiXiXi- 
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